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Abstract 



o. 

O , A layered antiferromagnetic superconductor in the mixed state may posses magnetic 

^ I domains created along the Josephson vortices. This may happen when an external 

^^' magnetic field is strong enough to flip over magnetic moments, lying in the phase 

c/2 . core of the Josephson vortex, from their ground state conflguration. The formation 

of the domain structure of the vortices modifies the surface energy barrier of the 
superconductor. During this process the entrance of the fiux is stopped and a newly 
created state exhibits perfect shielding. Such behavior should be visible as a plateau 
on the dependence of flux density as a function of the external magnetic held. The 
Q ' end of the plateau determines the critical held, which we call the second critical 

^/ field for fiux penetration. 



^ ■ Key words: High-Tc superconductivity, mixed state, magnetic superconductors 

0\ ; PACS: 74.60.Ec, 74.72.-h 
O 

in 

O ■ INTRODUCTION 

O' 

Among classical magnetic superconductors there are three groups of clus- 

rS ' ter compounds, REMogSg, REMoeScg, and RERh4B4( RE=rare-earth) which 

c^ ' have been the primary systems for study of the interplay between supercon- 



ductivity and long-range magnetic order [T||2] . Although good quality single 
crystals of these materials have been available and measured for a long time 
a very interesting phenomenon was recently discovered in DyMogSg due to 
very carefully conducted experiment [3]. This phenomenon predicted in |1] 
and named two-step fiux penetration was previously observed solely on [bet) 
ErRh4B4 0. The present work is inspired by this discovery and the hope that 
the same behavior could possibly be observed in some of the layered supercon- 
ducting structures. The specific feature caused by the long antiferromagnetic 
order in the mixed state of a superconductor is the possibility of creation of 
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the spin-flop (SF)(or inetamagnetic) domains along the vortices [6]. This is 
easy to understand by taking two sublattices antiferroniagnet as an example. 
An inflnitesimal magnetic fleld applied perpendicular to the easy axis makes 
the ground antiferromagnetic (AF) state unstable against the phase transfor- 
mation to the canted phase (SF). On the contrary, if the magnetic fleld is 
applied parallel to the easy axis the antiferromagnetic conflguration is stable 
up to the thermodynamic critical fleld Ht- When the fleld is further increased 
a canted phase develops in the system. Let us assume that in an antiferro- 
magnetic superconductor the lower critical fleld fulflls the relation Hd < \Ht 
and that the external fleld, Hd < H < \Ht , is applied parallel to the easy 
axis. Then the superconducting vortices appear in the ground antiferromag- 
netic state. When the fleld is increased above Hpi (see Figj3]) approximately 
equal to \Hj' the phase transition to the canted phase originates in the vortex 
core. The spatial distribution of the fleld around the vortex is a decreasing 
function of the distance from its center. Hence the magnetic fleld intensity in 
the neighborhood of the core is less then Ht- Therefore, the rest of the vortex 
remains in the antiferromagnetic conflguration. The radius of the SF domain 
grows as the fleld is increased. The above considerations apply to the classical 
superconducting Chevrel phases as well as to the high Tc superconductors, 
where antiferromagnetic order is produced by the regular lattice of RE ions 
occupying isolating layers. In this paper we consider the structure shown on 




Fig. 1. Schematic drawing of a piece of the layered superconductor. The shaded 
areas (n,n+l) represent superconducting layers. The bold arrows represent magnetic 
moments of RE ions lying in the isolating layers. The axes of the reference frame 
are shown. 



Fig. [Uthat we believe simulates a real structure of many antiferromagnetic lay- 
ered superconductors. A good candidate to show the above behavior should 
possess the isolating layers with the magnetic moments of RE ions running 
parallel and antiparallel to the direction (easy axis) lying in the ah plane. A 
typical example of such system is ErBa2Cu307. This compound has tetragonal 
unit cell with small orthorombic distortion in the ah plane. The Er ions form 
two sublattices antiferromagnetic structure of magnetic moments laying par- 
allel and antiparallel to the h direction [7] . Another example may be Er nickel 
boride-carbides Pl9p0] . The layered structure of RE nickel boride-carbides 
is reminiscent of that of the high- Tc oxide superconductors and consists of 



RE-carbon layers separated by Ni2B2 sheets. 



BASIC EQUATIONS 



We start description of our problem in terms of the Lawrence- Doniach en- 
ergy functional. In this approach a layered superconductor is described by the 
superconducting planes with the interlayer distance d, as shown on Fig. [TJ 
The antiferromagnetic subsystem consisting of RE ions is confined to the in- 
sulating layers. The magnetic moments are running parallel and antiparallel 
to the X-axis (easy axis). The Lawrence-Doniach functional is obtained from 
the standard Ginzburg-Landau energy by discretization of the kinetic energy 
in the z-direction. 
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The quantity %, e, m, denote Planck constant, charge of the electron and mass 
of the current carrier in the ab plane, whereas Ai denotes mass of the current 
carrier in the z-direction. The antiferromagnetic two sublattices subsystem 
with single ion anisotropy is described with the following energy density func- 
tional 
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where M„ = Mi^ -|- M2n is the sum of the magnetization vectors of the sublat- 
tices in the n-th insulating layer, M^^ is the component of the magnetization 
sublattice vector along the anisotropy axis in the n-th layer, J denotes the 
exchange constant between two sublattices, K is the single ion anisotropy 
constant, J\'y\ is the magnetic stiffness length, and Mq = |Mi„| = |M2„|. 
Since in the following we analyze the phenomena with characteristic length- 
scales much larger then the interatomic distance it is justified to omit the 
gradient term in J'm- The components of the total magnetization vector M 
have the following form in both sublattices: Miy = Mq sin 6i, Mix = Mq cos di, 
where 6i (canted spin angle) is the angle between the magnetization in the 
sublattice and the external magnetic field directed along the x-axis. The AF 
(^1 = 0, ^2 = tt) and SF phases (^i = —62 = 6) are in thermodynamic equilib- 
rium in an applied field equal to the thermodynamic critical field [TT] 



Ht = Mo[K{J - K)] 



1/2 



(3) 



The canted spin angle of the SF phase is then expressed as 

cose = —-—. (4) 

Finally we add the magnetic field energy to obtain the free energy of the entire 
system 

F = F, + j[fM + ^^']dV. (5) 

According to experiments the antiferromagnetic order is very weak affected 
by the presence of superconductivity, then it is reasonable to neglect the effect 
of superconductivity on the exchange interaction in F. Instead we introduce 
electromagnetic coupling between the magnetic and superconducting subsys- 
tem. This means that both order parameters \l/„ and M are coupled through 
the vector potential A 

B = rotA = /ioH + M, (6) 

js = rotH, (7) 

where B is the vector of the magnetic flux density (magnetic induction) and H 
is the vector of the thermodynamic magnetic fleld intensity. The functional ([5]) 
can be treated in the London approximation by assuming a constant modulus 
\l/„ within the planes and allowing only for phase ( (/?„) degree of freedom. The 
equilibrium conditions of the whole system can be obtained via minimization 

the Gibbs free energy functional G = F — (Blio)dV. Performing this task 

with respect to vector potential A and (pn provides us with the fundamental 
equations for currents and phases. 

E X^(^V(2)</^„ - A(2))5(;2 - nd) = /ioj(2) = rot(2)(B - M) (8) 

E (t^^^^^Xn+i,n)e{z - dn)e[d{n + l)-z]= /ioi. = rot,(B - M) (9) 



„ .271 X]d 

V(2)fV(2)V5n - — A(2)j = — UinXn+l,n-sinXn,n-lj (10) 

where 6{z — nd) is the Dirac delta function, Q{z — dn) Heaviside step func- 
tion, A denotes London penetration depth in the superconducting plane, Xj = 

I I 2ei /■{"+i)'^ 

XJM/m, r.j = dJM/m and Xn+i,n = V^n+i - fn + ^ A^dz is the 

n, Jnd 

gauge invariant phase difference. In the following we shall investigate the prob- 
lem of a single vortex line lying parallel to the Josephson coupled supercon- 
ducting layers, separated by the insulating antiferromagnetic layers. 




Fig. 2. Single Josephson vortex lying in the ab plane along the x-axis. The shaded 
area shows induced SF domain along the phase core. 

SINGLE JOSEPHSON VORTEX 



We assume that the vortex center, located in the central n = layer, is parallel 
to the X-axis. The relation between the magnetic field of the vortex and the 
gauge-invariant phase difference Xn+i,n{y) oi the superconducting wave func- 
tion across layers n and n + 1 can be found by integrating the vector potential 
given by equations (|8]) and ([9]) around a rectangular, semi-infinite contour C, 
as shown in Fig. [H This contour, parallel to the yz plane, is located apart 
from the central junction n = 0, where nonlinearities of the phase difference 
must be taken into account. The magnetic flux within this contour in given 
as. _ 

Hy) = d dy'B{y',z)= j Adl 

Jy JC 

Differentiating the result with respect to y one obtains 
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The above equation shows that in the Josephson vortex the screening currents 
vanish on a length scale Aj along y axis, and a scale A along z direction. On 
these scales the Josephson and Abrikosov vortices in an anisotropic supercon- 
ductor are roughly equivalent apart from small corrections in the current flow. 
Fig. [21 and the magnetic fleld pattern due to the layered structure. But in con- 
trast to the Abrikosov vortex, where the large current flow near the core leads 
to complete supression of the order parameter, the supression of the order 
parameter in the superconducting layers is only weak in the so called phase 
core of the Josephson vortex. Within the distances Tj along y, and d along 
z we have to take into account the nonlinearity and discretness of Eq. flTOj) . 
On these length scales the phase is changing rapidly and the current density 

jz reaches its maximum value jc 
London model fails. 
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In the region of the phase core 



To make the problem simpler we assume that the magnetization in the isolated 



vortex is constant across the SF phase domain. Thus, we can write 
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where pm is dimensionless radius of the SF domain in the cyhndrical reference 
frame x = x ; y = Xjp sin 6 ; z = XpsinO. Then the solution of Eq. (TTT]) for 
a single Josephson vortex is given by the modified Bessel functions Kq and Iq 



6sF = CiKo (p) + C2/0 (p) , for Pj < p< pn 
&AF = C'si^o (p) , for p > Pm, 



(13) 



{pj denotes the dimensionless phase coherence length ) with the following 
boundary conditions: 



&SF {pm) = PqHt + M = Bt 

bAF {Pm) = PqHt- 



(14) 



These conditions, together with the flux quantization condition, are used to 
calculate the arbitrary constants in Eq. flT^ . 
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Finally we write free energy of the isolated vortex 
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is the position of the vortex line, cxi denotes the surface 



of the phase core, and o"2 the surface of the SF domain respectively. The 
integrals in Eq. ( 1TB|1 performed as line integrals along the contours of the cross 



sections of the appropriate surfaces give ei - the hne tension of the vortex. 
The minimum of ei with respect to pm determines 

2 ^*^0 /IVN 



FREE ENERGY OF THE LATTICE 

Equation flTTl) . in the new coordinates, can be rewritten for the lattice of 
vortices in the following way: 



B + rotrotB = — — ^ 5{r — r 



AAj ^ 

where r^ specify the positions of the phase cores of the vortices. The solution 
of Eq. (fT8|) is then a superposition 

B(^) = J2'^rnir - rm) 

rn 

of the solutions Bm(^ — Tm) of isolated vortices at points rm- The free energy 
of the system can thus be written as 

AA r 

F = ^ f da(B X rotB) (19) 

The above symbolic surface integral is taken over the surfaces of the phase 
cores and surfaces of the SF domains. The energy of the Meissner state is 
chosen as zero of the energy scale. Again, when the surface integrals are re- 
placed by contour ones over appropriate cross sections we get line energy of 
the system. This, in turn, multiplied by vortex density n gives /-free energy 
density of the system. After some transformations one can derive the following 
formula 

f = ne,+ n<f)oHT{\nl3)-' Y. Mrm), (20) 



here the sum is over all vortices excluding the one in the origin, and r^ denotes 
the distance of a vortex from the origin. The lattice sum may now be replaced 
by integral in the yz-plane over a smoothed vortex density, excluding the area 
n~^ associated with the single flux line in the origin. The free energy density 
then reduces to 



/ = „.,+Bn^ ^i^ +B-i^,a>.>K^ (21) 
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here a = \ai\ denotes the length of the basal vector of the nonequilateral 
triangular unit cell, and 2\a2\ = aVl + tan^ a ( a is the angle between both 
vectors), tana = ./|^ [12]. To determine the equilibrium state it is necessary 
to minimize the Gibbs free energy density with respect to magnetic induction. 
The result yields an implicit equation for the constitutive relation B = B{H) 

H-"-^ = Br-^]{^^]+ii:^Ji;^\n{^=] (22) 




FLUX PENETRATION 



Consider semi- infinite specimen in the half space y > 0, the vortex and the 
external magnetic field running parallel to the surface in the x direction. The 
presence of a surface of the superconductor leads to a distortion of the field 
and current of any vortex located within a distance of the order of penetration 
depth from the surface. To fulfill the requirement that the currents cannot flow 
across the surface of the superconductor we need to introduce an image vortex, 
with vorticity opposite to the real one. Both vortices, direct and image, interact 
as real ones except that the interaction is attractive. In the low flux density 
regime, Clem [T3j has shown that there exist two regions: a vortex-free region of 
the width i/ff near the surface of the sample, and a constant flux density region 
for y > yff. Within the vortex- free area one can introduce the locally averaged 
magnetic field Bm which is a linear superposition of the Meissner screening 
field, the averaged direct vortices flux density exponentially decreasing towards 
the surface from its interior value B at y = yjf, and averaged image vortices 
flux density. In our problem the x component of this superposition can be 
approximated by 

Bm = B cosh ( ^^^~ A . (23) 

The boundary condition B^f{0) = ^qHq determines the thickness of the 
vortex-free region 

,„ = A,.cosh-'(^). (24) 

We assume that the test vortex line is lying within vortex free region at a 

y y 

point r = (-—,0), and its image at r = (— — ,0) outside the superconductor. 
Now the local field of the test vortex can be understood as a superposition of 



the following fields 



BsF = bsF (r) - bAF (2r) + xBaj {rjf - r) 
Baf = bAF (r) - bAF (2r) + xBaj {rjf - r) 

'Vff 
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where rjf = (-p-, 0), and x denotes the unit vector in the x direction. Having 

determined the local magnetic field we can write the Gibbs free energy of the 
test vortex line as 
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After some transformations p^jT^ one can obtain the Gibbs free energy per 
unit length Q 

G = Gi + G[ + gM, (27) 

where 
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Gi describes the interaction of the test vortex with its image, G'l is a correction 
term introduced by Clem |T3] , and Gm describes the interaction energy of the 
test vortex with the modified Meissner field. To find the conditions of the 
vortex entrance and exit, one has to solve a force balance equation for the test 
vortex, at the surface of the sample, and at the edge of the fiux- filled area, 
respectively. A calculation using Gi and Gm alone gives non vanishing force on 
the test vortex at r = Vff. However, the force should be zero there, because 
Gm is supposed to account for all the image vortices. To avoid double counting 
the image vortex one can subtract from the self-energy a contribution of the 



excess image fixed at r = —Vff. One can easily check that Q[ is neghgible at 
the surface of the sample and has no influence on the conditions of the flux 
entrance. When the flux starts to enter the sample, Hq = iJe„2(-B), 

Vff = yen = Xj cosh (^ j , (30) 

and the energy barrier is moved toward the surface within pm- Thus, one can 
derive from the force balance equation 



Di dbApip) 



2D2 dp 



= Bsmh(^]. (31) 
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The left hand side of the above equation gives ifen2(0) = -ffT/3(21n/3)~^. This 
field may be thought as the second critical field for flux penetration calculated 
in the single vortex approximation [H]. Combining Eqs. (!30|l and (l3T|l we 
finally obtain 

In the opposite case, when the flux exits the sample, the surface energy barrier 
tends to the edge of the flux-filled zone. Similar considerations as the above 
show that 

lJioHe.2{B) ~ B. (33) 

The measure of the height of the energy barrier against flux entrance is 

Ai7en(5) = |i/en2(5)-i7e,(5)|, 

and against flux exit 

Aif,,.(5) = |ife,(5)-iJex.2(5)|, 

where if eg is given by Eq. fl22|) . 



DISCUSSION OF THE RESULTS 



Let us make a short summary of the calculations and visualize the results on 
schematic magnetization curve shown in the Fig. [31 When the external field 
is not strong enough to create the SF domains inside vortices, than the mag- 
netization process of the sample being entirely in the AF phase is as follows. 
The vortices without magnetic structure start to enter the specimen at Hf.ni- 
When the field is increased up to the value Hpi, which is of the order of Ht, 
the SF domains are created. Now, the screening current must redistribute its 
flow in order to keep constant the flux carried by the vortex. This feature 
is easily seen from Eqs. (113] - [T5l) . The redistribution of the screening current 
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Fig. 3. Schematic drawing of the magnetization process, -ffeni denotes the first 
penetration field for vortices without magnetic structure, Hpi is the apphed field 
which originates SF transitions inside vortices, and Bpi is the corresponding flux 
density. Hen2 is the entrance field for the vortices possessing magnetic structure. 



changes the surface energy barrier preventing vortices from entering the sam- 
ple as expressed in Eq. (1281) . It means that the density of vortices n is kept 
constant. Consequently the averaged flux density in the sample B = mpo re- 
mains constant when the external field is increased. In Fig. [3] this feature is 
visible as a plateau on the B{H) curve, or alternatively as a second negative 
slope on the M{H) curve. The vortices start to penetrate the sample when the 
external field reaches the right edge of the plateau. We call this value, given 
by Eq. fl52]) . second critical field for flux penetration Hen2- 



To find the thermodynamic critical field Ht, and then to calculate Hen2{B) the 
following argumentation is proposed. At low fields, in the vicinity of the lower 
critical field ifci, the intensity of the field in the vortex core is 2Hci [I5]. When 
the external field is increased the field intensity in the vortex core increases 
because of the superposition of the fields of the surrounding vortices. The field 
intensity in the core must reach Ht in order to originate a transition to the 
SF phase. Thus, taking into account only the nearest neighbors we can write 
for the nonunilateral triangular lattice 



Ht = 2Hci + z 



^0 



nXXjfiQ 



Ka -\+2Kn — 



X 
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3A 

a; 



(34) 



here a corresponds to the value Bpi of the flux density for which the penetration 
process stops, see Fig. ([3]). From the relation Bpi = 2ipoJXj/{a'^\/3X) one can 
compute a, which in turn may be inserted back into Eq. ( 134|) . It is easy to 
estimate the saturation magnetization Mq taking into account the volume of 
the elementary cell. Then, Eqs. ([3]) and (jl]) can be used to calculate M in the 
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SF-phase domain 

M = 2Mocose= ° . (35) 



CONCLUSION 



The layered antiferromagnetic superconductor may reveal below T^r a very 
interesting behavior in the magnetic field applied parallel to the supercon- 
ducting planes. When the sample is in the virgin state, initially it magnetizes 
like ordinary type II superconductor. Upon the applied magnetic field of in- 
tensity equal to the critical field for fiux penetration the sample undergoes 
a transformation from the Meissner to the mixed state. Then, the magneti- 
zation may proceed in an unusual way. When the field is further increased a 
new state may appear in which vortices possesses the spin-fiop phases created 
around the cores. We have assumed that in this new state vortices undergo 
metamorphosis to the shape shown in Fig. |2l This state is characterized by 
the plateau on the magnetization curve, shown in the Fig. |3l It means that 
the magnetic fiux density inside the sample is unaffected by an increased ex- 
ternal field. This perfect shielding should occur until the applied field reaches 
certain value of intensity, we call it second critical field for fiux penetration. 
Then the vortices possessing magnetic structure enter into the sample. This 
phenomenon we named two-step fiux penetration. 
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